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The intrinsic alignment of galaxies constitutes the major astrophysical source of system- 
atic errors in surveys of weak gravitational lensing by the large-scale structure. We dis- 
cuss the principles, summarise the implementation, and highlight the performance of two 
model-independent methods that control intrinsic alignment signals in weak lensing data: 
the nulling technique which eliminates intrinsic alignments to ensure unbiased constraints 
on cosmology, and the boosting technique which extracts intrinsic alignments and hence 
allows one to further study this contribution. Making only use of the characteristic depen- 
dence on redshift of the signals, both approaches are robust, but reduce the statistical power 
due to the similar redshift scaling of intrinsic alignment and lensing signals. 



1 Cosmic shear and intrinsic alignments 

The weak gravitational lensing of distant galaxies by the large-scale structure, or cosmic shear in short, is one of the most 
powerful cosmological prob es of dark matter, dark energy, deviatio ns from general relativity, or the in itial conditions 
of structure formation (e.g. | Peacock et al 1 1200(4 lAlbrecht et al.ll2006h . Recent observational results (e.g. Benjamin et al 



_OH ... 
planned survey such as LSSlQ and Euclid^ feature cosmic shear as one of their primary probes. 

The gradually tightening statistical constraints on cosmology by upcoming cosmic shear surveys entail more and 
more stringent requirements on the control of systematic errors. While instrumental effects, which can e.g. affect the 
measurement of gal axy shapes, ar e eliminated via specific instrument designs, survey strategies, or steps in the reduction 



2007; Schrabback et al. 2010) demonstrate the increasing maturity and competitiveness of cosmic shear, while ambitious 



pipeline (see e.g. Bri dle et al . 2010), astrophysical sources of systematics have to be dealt with at the data level. The major 



astrophysical systematic in cosmic shear surveys is generated by the intrinsic alignment of galaxy shapes, which we will 
consider in detail in the following. 

At the two-point level cosmic shear measures, such as the correlation functions or the power spectrum, are based 
on t he correlator of the (complex ) measured ellipticity e of pairs of galaxies (for an overview on weak lensing theory 
see iBartelmann & Schneiderll200ll) . In the limit of very weak lensing this ellipticity is given by the sum of the intrinsic 
ellipticity of a galaxy e s and the gravitational shear y. Consequently, a two-point correlator of galaxy samples i and j has 
contributions by four terms, 

H) = (y^)-(^r)-(y^X^y;)- d) 

GG II GI 

where the first term on the right-hand side is the desired cosmic shear signal (GG henceforth). The second term denotes 
correlations between the intrinsic ellipticities of galaxies (II henceforth), and the third term comprises correlations between 
the intrinsic ellipticity of one galaxy and the gravitational shear acting on another (GI henceforth). In standard cosmic 
shear analysis it is assumed that the intrinsic shape of a galaxy is not correlated with either the intrinsic shape or the 
gravitational shear acting on another galaxy. However, galaxies do indeed show intrinsic alignments, which thus induces 
systematic errors. 

A straightforward solution to this issue would be the modelling of the II and GI signals, possibly including a set 
of nuisance parameters in the cosmic shear likelihood analysis, but this approach is only safe if the functional form 
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of the systematic is well known (e.g. iKitching et al. 2009). Unfortunately, the understanding of intrinsic alignments is 
currently still at a crude level because it involves the intricacies of galaxy formation and evolution in the local large- 
scale environment, including the complicat ions by baryonic physic s. Analytic modelling of intrinsic alignments is only 
applicable on relatively large scales (see ISchneider & Bridle! 2010l for rec ent developments), and simulations are eithe r 
limited by simplistic models of placing galaxies into dark matter haloes (e.g. Hevmans et al.l2006t Semboloni et al.l2 008). 
or by the small volume if baryons are included (e.g. lBett et al ] |2010l:lHahnet alll2010h . 

Observations of intrinsic alignments, using shal low surveys in which cosmic shear is subdominant or cross-correlations 
between galaxy num ber densities and shapes (e.g. iMandelbaum et al. 20061 : Hirata et al. 2007 : Mandelbaum et al. 2010; 



Joachimi et al. 2010h . have revealed that this systematic depends on galaxy colours (and hence type) as well as galaxy 



luminosity (and hence mass). However, constraints on specific intrinsic alignment models are still weak and inherit sys- 
tematic uncertainties from e.g. galaxy luminosity functions or ^-corrections. Furthermore, these observations are only 
possible for galaxy samples which are either luminous or at low redshift, and therefore can have spectroscopic or high- 
quality photometric redshift measurements, so that a substantial amount of extrapolation is needed to predict intrinsic 
alignment signals from the typically faint and high-redshift galaxies in cosmic shear surveys. 

Hence, model-independent methods that remove the contamination by intrinsic alignments from cosmic shear surveys 
are desirable. Besides, it would be convenient to be able to extract the intrinsic alignment signal from the cosmic shear 
survey itself in order to at least provide a cross-check for the extrapolation of intrinsic alignment models from brighter 
and low-redshift samples. These two goals are met by the nulling and bo osting techniques which will be presented in the 
following, based on the work by Joachimi & Schneidei (2008, 20091 2010s, Paper I, II, III hereafter). 



2 The principle of nulling and boosting 

As cosmic shear two-point statistics are based on the correlator (Q~|), all of them are affected by additive contributions from 
II and GI correlations, so that for the observed angular power spectrum^ one obtains 

c lbsW = c ggW + c rf W + c giV) • (2) 

where the superscripts in par entheses indicate that we consider tomographic cosmic shear data with additional photometric 
redshift information (see e.g. Schrabback et aljfeoiol) . correlating photometric redshift bin i with bin j. All power spectra 



on the right-hand side of (0 are given in terms of their underlying three-dimensional pow er spectra by the following 



Limber equations (e.g. lHirata & Seliakll20o4lBridle & King||2007l:|joachimi & Bridlef zOlO). 
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Here, P$ denotes the matter power spectrum, Pn the intrinsic ellipticity power spectrum, and Pgi the cross-power spectrum 
between the matter field and the distribution of intrinsic galaxy shapes. The latter two have to be specified by an intrinsic 
alignment model, but are generally unknown. Note that we assume a spatially flat universe throughout. 

The integrals in (0) to © run over the comoving distance^ out to the comoving distance horizon Xhoi- The kernels 
are composed of the probability distribution of comoving distances p^'fx) m photometric redshift bin i and the lensing 
weight 
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where the integral corresponds to the typical scaling of shear signals with the average over the source galaxy distribution 
i of the ratio of distance between lens and source over the distance between observer and source. The distributions p^(x) 
are usually chosen to be narrow in cosmic shear tomography, whereas the lensing weight (O is broad. These two quantities 
determine the scaling of the power spectra (O to (0 with redshift, which holds in particular for the II and GI signals, if 
P^ix) i s compact enough that the evolution of the three-dimensional intrinsic alignment power spectra over the integration 
range is negligible. 



3 Since the methods which are going to be discussed do not depend on angular scales, the following considerations hold for any cosmic shear 
two-point statistic. We use power spectra throughout because they are the easiest to handle computationally. 
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Figure 1: Redshift dependence of the GG, GI, 
and II signals. Shown is the amplitude of the 
power spectra C^' =1 '^(€) at i » 200 as a func- 
tion of the median redshift of galaxy sample j. 
The power spectra are computed for a ACDM 
cosmology and a Euclid-like redshift distribution 
divided into 20 bins containing an equal number 
of galaxies, and using a photometric redshift dis- 
persion of 0.05(1 + z). The GI and II terms are 
obtained from the non-linear version of the lin- 
ear alignment model. The GG signal is shown as 
black circles, the GI signal as dark grey triangles, 
and the II signal as light grey diamonds. 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

z 



As a first step, we illustrate the redshift dependence of the GG, GI, and II signals. To compute the matter power 
spectrum, we assume a flat ACDM universe wit h parameters Q m = 0.25, <r 8 = 0.8, n s = 1, h = 0.7, and Qb = 05, using 
the transfer function by Eise nstein & Hu (1998) and the non-linear corrections provided by Smith et al. ( 20031). The two 



intrinsic alignment power spectra are determined via the non-linear version of the linear alignment mod el (ICatelan et al 



200ll:lHirata & SeliakH 20Q4: Bridle & King 2007), with the normalisation given in lBridle & King! d2007l) . We emphasise 
that an intrinsic alignment model is only introduced for illustrational purposes and to assess the performance of the 
methods we propose, while the implementation of these techniques does not depend on any knowledge about intrinsic 
alignments other than the well-known form of and (|5]Q. Here and in the following we will assume a Euclid-like survey 
covering 2 0000 deg 2 with a m edian redshift of 0.9 and a photometric redshift scatter of cr ph (l + z) with <x ph = 0.05 (for 
details see lLaureiis et al l2009l) . 

In Fig.[TJwe have plotted the angular power spectra for i = 1 and I « 200 fixed as a function of the median 

redshift of the photometric redshift bin j. We have split the survey into N z = 20 photometric redshift bins such that each 
bin contains the same number of galaxies. It is evident that the II signal has a distinctive dependence on redshift, quickly 
decreasing in amplitude if the distance between bins i and j increases. This behaviour is related to the narrow kernel in (0]) 
and originates from the fact that galaxies have to be physically close in order to be aligned by the sam e matter structures. 
Thus the II signal is relatively straightforward to remove or iso late, also in a model-independent way ( King & Schneidej 
20021 120031: iHevmans & Heavensll2003l:lTakada & Whitdl2004l) . In contrast, the GI and GG signals display a very similar 
scaling with redshift as both have lensing contributions. In the following we will therefore concentrate on the GI signal, 
discussing the treatment of II correlations mostly in a qualitative manner. 

To construct methods that make use of the characteristic dependence on redshift of the GI and GG signal displayed 
in Fig.Q] we first consider the limit of very narrow redshift bins, i.e. p®(x) ~ 8d(x ~Xi)> where 5u is the Dirac delta- 
distribution, and where \i is a comoving distance on which p^'Xx) is centred. The lensing weight (0 can then be written 
as 



q (i) (x) -» q(xi,x) = 
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With this equation at hand, the angular GG and GI power spectra turn into 



Ccc(Xi,XjJ) 



- f 

Jo 



(xixj) qixi,x)q(xpx) „ (e 



dx 



xr 



q{XhXi) Xi 2 Pei [~>Xij + QOCuXj) Xj 2 Psi \~.<Xj 



(7) 



(8) 



(9) 



where we have made the dependence on^, and^ ; explicit in the arguments. Note that the first term in (0 only contributes 
if Xi < Xj' an d v i ce versa for the second term. Note furthermore that in this approximation the dependence of Cgi on the 
background bin (e.g. bin j in the first term) is only through the lensing weight q. 



Recent observational results iMandelbaum et al. 2010; Joachimi et al. 2010) suggest that the model we use over-predicts the intrinsic alignment 
signal of typical galaxies in deep cosmic shear su rveys. In addition, the implemented model has been shown in the meantime to contain a mathematical 
error (Hirata 2010, in prep.; Joachimi et al. 2010). Both issues do not have a significant impact on the interpretation of our results. 
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We now define transformed power spectra as simple linear combination of the original angular power spectra, summing 
over the index of the background photometric redshift bin, 



= 2 rgi^o, (io) 

j=jmm 

for every foreground bin i, where the sum starts at the bin index j m i n . We are going to determine the weights 7*® such 
that either the GI signal is eliminated ('nulled'), or the GG signal is suppressed (and GI correspondingly 'boosted'). 

Retaining the approximation of narrow redshift bins, one obtains by inserting (0 into ( fTUb that the GI signal is removed 
for / m j n > i if the weights fulfil the condition 
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Since auto-correlations with i = j are likely to be contaminated by an II term, we discard them completely and set 
7min = i + 1 . Then N z - i weights enter a linear combination (TT~0b which only has to fulfil the condition (fTTT i. so that several 
such linear combinations can be constructed, indexed by the subscript [m]. 

To obtain the central condition for boosting, we insert ([H) into (TlOb instead, finding that the GG signal is eliminated if 
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is satisfied. Here, we defined the function 

i=k V 

in a discretised form which can be used in ( TTZb if the integral is approximated by its Riemannian sum. To arrive at the 
second equality in ( fT2b . we inserted (0 and assumed that the redshift evolution of the matter power spectrum is given 
by Pg(k,z) oc (1 + zY 2 which holds in the matter-dominated era. Since the matter power spectrum is already known to 
fairly high precision, one could use the first term in (TTZb directly, but we choose to work with the second term to ensure 
a maximum of model independence. Since the weights 7^ define G (l) for a set of arguments Xk, there is no additional 
freedom (except for the arbitrary normalisation of the weights) to construct further linear combinations, contrary to the 
case of nulling. Hence we will use a single set of weights to extract the GI signal in redshift bin i. 

Note that the weights depend on comoving distances via ( fTTT > and ( TOl l. while the observed signals Cr'P(i) are binned 
in terms of redshift. Therefore the distance-redshift relation and consequently cosmological parameters like Q m , or the 
dark energy equation-of-state parameters, enter the weights 7^, but this dependence is weak, and even grossly incorrect 
a-priori assumptions for the values of these parameters are readily corrected for in an iterative scheme (see Paper II). 



3 Nulling intrinsic alignments 

The additional freedom in the determination of the nulling weights (fTTb can be exploited to maximise the information on 
cosmology contained in the transformed measures (TTOb . As investigated in Paper I, this can be achieved by calculating the 
weights such that a functional of the Fisher matrix like the trace or determinant is maximised. The ch oice of fu n ctiona l 
or Fisher matrix element does not have a strong effect on the resulting form of the weights, see e.g. Shi et al. (2010). 



Higher-order modes of the transformed power spectra nj^j(^) with m > 1 are then obtained by requiring that the corre- 
sponding weights be orthogonal with respect to the weights of all lower-order modes. We showed in Paper I that nearly all 
cosmological information is already contained in the first order m — 1. In Fig.|2]an exemplary set of weights T^l for three 
different bins i has been plotted. Due to the nulling condition (fTTb . these weights necessarily have a zero crossing. They 
are largest if foreground bin and background bin are close in redshift where the GI signal is small, and have the smallest 
absolute values at high redshift were the GI term is large relative to the lensing contribution. 

Alternatively, one can drop the condition on maximising cosmological information and simply compute sets of weights 
that are mutually orthogonal. For N z - i power spectra entering (TTOb for a given i there are just N z - i such sets of weights 
where one, however, cannot fulfil the condition (fTTb anymore, so that the corresponding transformed power spectrum is 
discarded in the further analysis. Imposing in addition unit normalisation on these sets of weights, this approach allows 
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Figure 2: Example set of nulling weights which 
additionally have been optimised to retain a 
maximum of cosmological information. The 
weights have been obtained for an equidistant 
redshift binning with N z = 20 between z = 0.2 
and z = 4. Weights are shown for the foreground 
bins ; = 1 (in black), i = 3 (in red), and i = 5 (in 
blue), with the lower redshift boundary of these 
foreground bins given in the legend. 
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for the convenient interpretation that nulling is equivalent to a rotation of the cosmic shear data vector, with a subsequent 
truncation of certain contaminated elements in the rotated data vector. This ansatz was investigated in Paper II and is 
computationally preferable because it does not require Fisher matrices to determine the nulling weights, but only (fTTt and 
a standard orthogonalisation procedure. Note that in this case no particular order is assigned to the modes m. 

In Fig. [3] left panel, the effect of the nulling transformation on power spectra is demonstrated, assuming a Euclid-like 
survey as detailed in the foregoing section, but for N z — 10 and a lower photometric redshift scatter cr p h = 0.020 While 
for the original power spectra C^(^) the absolute value of the GI signal attains about 10 % (for auto-correlations, i = J) 
to 50 % (for a high redshift of the background bin, j = N z ), the GI term is suppressed by roughly two orders of magnitude 
with respect to GG in the corresponding nulled power spectrum. Depending on the mode m, the residual GI signal can 
feature one (as in the depicted case) or more zero crossings and hence acts more like a noise contribution (see Paper II). 

The GI signal would be perfectly eliminated if the redshift bins were infinitesimally narrow. In a finite size photometric 
redshift bin, the redshift zu corresponding to the comoving distance \i where the intrinsic alignment contribution is 
removed, has to placed such that the GI term is still largely downweighted. Simply choosing the mean of the photometric 
redshift bin boundaries has proven to yield a robust and efficient nulling method, see Paper II for details. Furthermore, 
if the distributions p®(x) have finite width, the second term in (f5J) also yields a contribution, even for i < j, because due 
to the photometric redshift scatter a galaxy from the background bin j may actually be located in front of galaxies from 
bin /. This additional GI contamination is not accounted for by nulling but is limited to photometric redshift bins that are 
adjacent in redshift. We suppress this term by downweighting contributions of power spectra C^(€) to ( TTOb . the stronger 
the smaller j - i, as detailed in Paper II. Note that this measure conveniently removes II correlations in the presence of 
photometric redshift scatter as well. 

Switching back to a photometric redshift scatter cr ph = 0.05 which should be achievable for most galaxies in future 
cosmic shear surveys, we predict the statistical constraints and the bias due to intrinsic alignments on a Euclid-like 
survey divided in to N z = 20 photometric redshift bins. To this end, we perform a Fisher matrix analysis on the set of 
cosmological parameters {Q m , cr 8 , h, n s , Q.b, Wo, w a ], using Gaussian covariances for the power spectra that account for 
both cosmic variance and shape noise. The resulting 2<x confidence contours, for the case where intrinsic alignments 
are ignored, and after applying the nulling technique, are shown in Fig. [4] The strong intrinsic alignment contamination 
we have assumed in this study is reliably reduced to insignificance by nulling, including the one by the II term which 
is downweighted by the measures described above. The price to pay is a considerable loss in constraining power which 
amounts to an increase of the lcr marginalised errors by factors of two to three on the cosmological parameters. Since 
nulling relies on redshift information to suppress intrinsic alignments, the dark energy equation-of-state parameters, which 
are primarily probed via the redshift evolution of structure growth and spacetime geometry, suffer most. It was proven in 
Paper II that this information loss is inevitable when summing over redshift bins as in ( TTOb . even if precise redshifts are 
available, and is a consequence of the similar redshift dependencies of the GI and GG signals, see Fig.Q] 

Since the nulling technique only makes use of the characteristic redshift dependence of the lensing and intrinsic 
alignment signals, its performance strongly depends on the photometric redshift accuracy. It was demonstrated in Paper 
II that nulling places stringent requirements on the scatter <x ph and the fraction of catastrophic outliers in the relation 
between photometric and spectroscopic redshifts (although the nulling variant employed in this section is robust up to at 
least <T p h = 0.1 and 6 % catastrophics, discarding only slightly more cosmological information if the photometric redshift 

5 Note that t he analysis in th i s secti on is based on a slightly different cosmology with erg = 0.9, using the non-linear corrections to the matter power 
spectrum by Peacock & Dodds 1 1996). These differences should have a very minor influence on the compatibility of our results. 



5 



j=3 







(,(, - 




GI 


■ 1 | 1 H-l- 

j=6 


1 — i 




GG - 




GI 


, — i 1 — 

7=10 


^-++-H 1 1 : 




GG - 


— — i — i 1 — i—i- 


GI 

1 1 : 


nulled 

_ 




GG - 't 




GI 



0.001 
le-04 



cr le-05 
O 

le-06 



le-07 
0.001 

le-04 



le-05 

le-06 
0.001 

le-04 

le-05 



le-06 
le-04 
le-05 
C le-06 
le-07 
le-08 
le-09 



GG 
GI 



mid 



GG 
GI 



far 



GG - 
GI 



boost 



GG 
GI 



10 



100 



1000 



10000 



10 



100 



1000 



10000 



Figure 3: Left panel: Effect of nulling on GG and GI power spectra for a Euclid-like survey divided into N z — 10 bins, 
each containing the same number of galaxies, and with o" p h = 0.03. Shown are the original GG (solid lines) and GI 
(dotted lines) power spectra for foreground bin i = 3 and background bin 7=3 (top panel), j = 6 (upper centre panel), 
and j — 10 = N z (lower centre panel). The bottom panel displays the nulled power spectra for i = 3 and m = 1. Right 
panel: Effect of boosting on GG and GI power spectra, also for cr ph = 0.03, and obtained for the default binning used. The 
coding of the curves is the same as above. The upper three panels show power spectra for a foreground bin at z, = 0.76 
and different background bins j, i.e. auto-correlations (j = i, 'auto'), cross-correlations with a bin at intermediate redshift 
(j = (i + N z )/2, 'mid'), and cross-correlations with the most distant bin (j = N z , 'far'). In the bottom panel the boosted 
GG and GI signals are plotted. Note that absolute values of the power spectra are shown throughout. 



quality deteriorates), as well as the uncertainty in the mean of redshift di stributions. Howev er, these requirements are met 
by - or actually drive - the goals of planned cosmic shear surveys (e.g. lLaureiis et al ] l2009h . 
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Figure 4: Parameter con- 
straints before and after nulling. 
Shown are the two-dimensional 
marginalised 2<x-errors for the 
full data set (if intrinsic align- 
ments are ignored) as blue 
contours and for the nulled data 
set as orange contours. The fidu- 
cial parameter values are marked 
by the crosses in each panel. 
The Euclid-like survey has been 
divided into N z = 20 photometric 
redshift bins with <x ph = 0.05. 
The parameters {h, Cl\,,w a } have 
been marginalised over. We 
employed our default intrinsic 
alignment model. 



4 Boosting intrinsic alignments 

Since in the case of intrinsic alignment boosting the weights in (TTOb are used to define the function ( fT3l ) over the integration 
range of (TT2l) . a dense binning in redshift and hence a large number of weights 7^® . is appropriate. As one does not gain 
information if the redshift distribution is sliced into bins which are much finer than the typical redshift scatter, the quality 
of redshifts drives the maximum useful number of bins and is thus again paramount for the performance of this technique. 
We consider three galaxy samples, one with a standard photometric redshift quality of <x p h = 0.05 ('P2'), one with high- 
quality photometric redshifts cr p h = 0.03 ('PI'), and one with spectroscopic redshifts assuming <x ph = ('S'), using all of 
them in combination with realistic number densities as expected for a Euclid-like survey (see Paper III for details). 

We use photometric redshift bins with a separation of 0.01(1 + z) and 0.02(1 + z), respectively, so that one can to good 
approximation replace the set of weights . with j — 0, .. , N z by a smooth weight function termed B®(x). In Paper III it 

was shown that this weight function is readily determined from the function G®ix) in ( fT3l via B^(x) = ^ , d 2 G (,) (y)/d^ 2 . 
The function G w (y) in turn is chosen such that ( TTZb is fulfilled and at the same time the transformed GI signal remains 
strong. Among many options this can be achieved by using a parametric form 

G (i) (x) = N expf- ^'f^ 2 ) (x-b), (14) 



where N, <x, b, and Xm are free parameters. These are fixed by requiring <TT~2t . a maximum at z, = z(xd within bin i, 
and a fixed (arbitrary) normalisation. An example weight function for the sample PI and z; = 0.76 is presented in Fig. [5] 
The Gaussian in dT4b dominates the peak around z, while the term^ - b allows for the negative minimum in G (,) (x) at 
low redshift that is necessary to fulfil d!21 l. Since the corresponding weight function B^(x) varies considerably over the 
redshift range, a dense sampling by photometric redshift bins (indicated by the crosses) is required. 

For the same set of parameters and survey specifications we have plotted original and boosted power spectra in Fig. [3] 
right panel. While the power spectra C ( ^\f) have a similar level of contamination as for the case shown in the left panel 
of the figure, the transformed power spectrum is now clearly dominated by the GI signal, the lensing contribution being 
suppressed by more than an order of magnitude with respect to the GI term on all angular scales. 

To assess the signal-to-noise of these boosting-transformed power spectra, we again resort to a Fisher matrix analysis, 
this time constraining parameters of a toy intrinsic alignment model and treating the residual contributions by cosmic 
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Figure 5: Functions G w (y(z)) and B (c, (x(z)) for 
Zi = 0.76, as indicated by the vertical grey line. 
In addition the sampling points corresponding to 
the median redshifts of the bins are shown for 
B^ixiz))- The width of the Gaussian has been 
chosen to cr z = crdz/dx — 0.085. Note that the 
normalisation of the functions is in principle ar- 
bitrary. 




Figure 6: Constraints on the free parameters of the GI model. 
Shown are the 1 cr confidence contours for three different survey 
models. The solid ellipse with its centre indicated by a filled circle 
corresponds to the spectroscopic survey S, the dashed ellipse and 
open circle to survey PI, and the dotted ellipse with diamond to 
survey P2. Note that the centres of the contours are offset due to 
the bias by the residual GG signal. The cross marks the fiducial 
values of A — 1 and y = 0. 
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shear as the systematic that causes a bias on these parameters. We set up a two-parameter model of the form P l ™ dd (k, z) 
A P S i(k, z) [(1 + z)/(l + Zo)V with the two free parameters A and y, where zo is a pivot redshift and Pgi(k, z) is given by 
our default intrinsic alignment model described in Sect.|2] In Fig.|6]the lcr confidence regions for the intrinsic alignment 
model parameters are displayed. For all three galaxy samples that we consider the bias due to the residual GG signal 
is clearly subdominant to the statistical errors. Constraints by the spectroscopic sample S are very weak due to the low 
number density of galaxies for which spectra are available; those for the two photometric redshif t samples are similar and 
comparable to curre nt constraints from galaxy number density-shape cross-correlations (e.g. Mandelba um et al. 2010; 
Ijoachimi et al.|[2oTob . 

However, current surveys are much smaller than the one assumed for our Fisher analysis, so that the parameter con- 
straints appear surprisingly weak. As we could demonstrate quantitatively in Paper III, the low signal-to-noise in intrinsic 
alignment boosting is directly related to the severe loss of information when applying the nulling technique, both governed 
by the difficulty of separating the GI and GG signals only by means of redshift information. Note that we did not include 
the II signal into the calculations of this section. Since the weight function B (,, (x) peaks close to Zi, see Fig. [5] the II term 
would also be boosted with respect to the lensing signal. If one preferred to study the II and GI signals individually, the 
boosting transformation could be preceded by the extraction of physically close pairs of galaxies from the sample, where 
the efficiency of this additional step again depends critically on the quality of redshift information. 



5 Conclusions 

We demonstrated that by means of purely geometrical transformations of cosmic shear two-point statistics it is possible 
to both remove and extract the intrinsic alignment signal from lensing data. Both techniques robustly isolate the desired 
signals in a model-independent way and work for survey requirements compliant with upcoming large-area tomographic 
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cosmic shear surveys. Since both methods make use of the characteristic dependence of the intrinsic alignment and lensing 
signals on redshift, the demands on (photometric) redshift quality are stringent. The performance of these approaches 
is fundamentally limited by the similarity of the GG and GI signals as a function of redshift, necessarily implying a 
considerable loss of cosmological information in the case of nulling and a low signal-to-noise in the boosting-transformed 
intrinsic alignment signals. 

This latter finding suggests that the nulling and boosting techniques are unlikely to become the methods of choice 
to clean the lensing signal from intrinsic alignments or characterise intrinsic alignment signals in the galaxy sample 
for future cosmic shear surveys with ambitious goals on cosmological parameter constraints. However, due to their 
model independence and robustness, these methods can nonetheless serve as valuable cross-checks for intrinsic alignment 
removal methods that rely on stronger assumptions about the contamination, as well as for extrapolations of intrinsic 
alignment models to the redshifts and luminosities relevant for cosmic shear. 

If one wishes to keep the loss of statistical power due to intrinsic alignment control in cosmic shear analyses to a 
minimum, one can either mode l the II and GI signals , incorporating a-priori information about the dependence on redshift 
and/or angular scales (see e.g. Bridle & King 2007), or make use of external data to recover parameter constraints, in 



particular from galaxy number density correlat ions and cross-correlations between gala xy number density and shapes 
which come for free with cosmic shear surveys dBernsteinll2009l:ljoachirm & Bridldl201oh . 



Using simple combinations of data like in (fTUb has been identified before as an efficient, model-independent approach 
to eliminate potential sources of systematic errors, e.g . to remove small-scale information from the large-scale structure 



■ I i 

affected by baryonic physics (Huterer & White 2005) or to remove foreground contamination from cosmic microwave 
background multi-frequency data feennett et alJl2003h . The potential benefits from the se techniques for cosmic shear 
and related cosmological probes have not been exhausted with the present work, e.g. Shi et al.l ( 2010) have applied 
nulling to cosmic shea r three-point statistics for which the contamination by intrinsic alignments is potentially severe 
( Sem boloni et al . 2008), )ut at present completely unknown. Extensions that allow for the inclusion of prior information 



on the intrinsic alignment signal, or that make use of nulling to extract the lensing magnification signal from galaxy 
number density-correlations are currently under investigation. 
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